Fractal Weyl law for Linux Kernel Architecture 
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We study the properties of spectrum and eigenstates of the Google matrix of a directed network 
formed by the procedure calls in the Linux Kernel. Our results obtained for various versions of the 
Linux Kernel show that the spectrunr is characterized by the fractal Weyl law established recently 
for systems of quantum chaotic scattering and the Perron-Frobenius operators of dynamical maps. 
The fractal Weyl exponent is found to be v w 0.65 that corresponds to the fractal dimension of 
the network d « 1.3. An independent computation of the fractal dimension by the cluster growing 
method, generalized for directed networks, gives a close value d ~ 1.4. The eigenmodes of the 
Google matrix of Linux Kernel are localized on certain principal nodes. We argue that the fractal 
Weyl law should be generic for directed networks with the fractal dimension d < 2. 

PACS numbers: 89. 75. He, 05.45. Mt, 07.05.Bx, 05.45.Df 



I. INTRODUCTION 



The celebrated Weyl law [T[ determines the number of 
quantum states in a given energy interval of conservative 
Hamiltonian systems in a phase space dimension d as a 
function of the Planck constant h 



N E oc ff 



d/2 



(1) 



Recently this law has been generalized to the fractal Weyl 
law applicable to open quantum systems with complex 
energy spectrum of nonunitary quantum operators [2- 
7]. In this case the exponent v is linked to the frac- 
tal dimension of invariant sets of classical non-escaping 
orbits. Such type of operators appear in physical prob- 
lems of quantum chaotic scattering, open quantum maps, 
metastable quantum states in molecules and nuclei. 

In addition to open quantum systems it has been 
shown that the fractal Weyl law is also valid for the 
Perron-Frobenius operators of classical dynamical sys- 
tems [S| with absorption or dissipation [Qj]. A finite size 
matrix approximant S of such operators is efficiently con- 
structed using the Ulam method [10| which divides the 
whole phase space on N cells of size e. The eigenval- 
ues Xi of the Ulam matrix approximant S of the Perron- 
Frobenius operator are located in a complex plane of a 
circle with |A| < 1. For two-dimensional (2D) maps we 
have the matrix size N ex 1/e 2 so that a cell area e 2 ef- 
fectively plays a role of h and N\ ex N v where N\ is a 
number of eigenvalues Xi with |A| < |Aj| < 1. For dynam- 
ical strange attractors and 2D maps with absorption we 
have v = d/2 where d is a fractal dimension of a dynam- 
ical invariant set, e.g. a fractal dimension of a strange 
attractor [9(. For v < 1 the majority of eigenvalues have 
|A| -> for N ->• co. 

In fact the Ulam method applied to ID and 2D maps 
naturally generates the Ulam networks [ll] which prop- 
erties can be rather similar to the properties of scale-free 
directed networks such as the World Wide Web (WWW). 
Such type of networks naturally appear in various fields 



of science (see e.g. |12l - ll8l |). Therefore we can expect 
that the fractal Weyl law can appear also for directed 
networks similar to those of the WWW. The properties 
of such networks are well characterized by the Google 
matrix G [19| constructed on the basis of directed links 
between nodes: 



Gij = aSij + (l-a)/N 



(2) 



where the matrix S is obtained by normalizing to unity all 
columns of the adjacency matrix, and replacing columns 
with zero elements by 1/N, N being the network size 
[Lq . In the WWW context, the damping parameter a 
describes the probability to jump to any node for a ran- 
dom surfer. For < a < 1 the only eigenvector with 
A = 1 is the PageRank vector which has non-negative 
components and plays an important role in ranking of 
the WWW sites [H, GJ, d • The matrix G belongs to 
the class of Perron-Frobenius operators [lg|. The spec- 
trum of G for the WWW university networks was stud- 
ied recently in [20j and it was shown that a significant 
fraction of eigenvalues is concentrated at |A| — > that is 
a characteristic property of spectra with the Weyl law. 
However, no direct evidence for the fractal Weyl law for 
the WWW was found there for the reasons we explain 
below. 

It is interesting to note that many complex networks 
are self-similar and have fractal properties (see e.g. 
[2lL \22\ and Refs. therein). Thus it was shown that 
the fractal dimension of networks of WWW and actors is 
d = 4.1, 6.3 respectively [21(. It is argued that the origin 
of fractal architecture of network is linked to a strong ef- 
fective repulsion between the most connected nodes [22J • 
Such a generic fractal structure of complex networks also 
indicates that under certain conditions the fractal Weyl 
law can appear in the spectrum of the Google matrix 
of such networks. However, it is also important to note 
that usually the fractal properties are investigated for 
undirected networks or by converting directed networks 
into undirected ones (see e.g. [2l|, |22J). Here, on a con- 
crete example of the Linux Kernel network we show that 




FIG. 1: (Color online) Distribution of eigenvalues A in the 
complex plane for the Google matrix G of the Linux Kernel 
versions 2.0.40 (top left panel), 2.4.37.6 (top right panel) and 
2.6.32 (left bottom panel) with N = 14079, N = 85756 and 
N — 285509 respectively. Bottom right panel shows eigen- 
value distribution for the Google matrix G* with inverted 
link directions in the PCN of the Linux Kernel version 2.6.32 
with TV = 285509. Solid lines represent the unit circle and the 
lowest limit of computed eigenvalues in top right and bottom 
panels. 



the directionality of the network places a crucial role so 
that the directed network has different fractal dimension 
compared to the converted undirected network with the 
same nodes. The importance of directionality is very 
clear from the spectrum of the Google matrix: the spec- 
trum is in a complex plane for a directed network while 
for an undirected network the spectrum is on a real line. 
The paper is composed as following: the specrum of 
the Google matrix is analyzed in Section II, the fractal 
properties of the Linux Kernel network are considered in 
Section III, the properties of Google matrix eigenstates 
are studied in Section IV and discussion of the results is 
given in Section V. 



II. SPECTRUM OF GOOGLE MATRIX 

To check the validity of the fractal Weyl law for a cer- 
tain type of scale-free networks we focus on the Proce- 
dure Call Network (PCN) of the Linux Kernel software 
introduced and studied recently in [23J. In the PCN the 
directed links are formed by the calls between procedures 
and the Google matrix @ is constructed in the same way 
as for the WWW. The results of [23| show that only about 
1% of eigenvalues have |A| > 0.1 that is by a factor 50 
smaller compared to the university networks of similar 
size studied in [20] - This gives a significant indication 
on appearance of the fractal Weyl law in the PCN. To 




FIG. 2: (Color online) Dependence of the integrated density 
of states W on the relaxation rate 7 = —2 In (A) for the Google 
matrix of different versions of Linux Kernels from VI. with 
N = 14079 to 2.6.32 with N = 285509 as it is shown on 
the legend. The highly degenerate values of 7 are marked by 
arrows at 7 m = —2 In (a/m) for m = 1, 2, 3, 4, 5, 6. 



study the spectral properties of G we use the data for the 
Google matrix at a = 0.85 obtained in [2j]. Our results 
are not sensitive to a which affects the spectrum of G in 
a rather simple way [IJ, [l6| . 
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FIG. 3: (Color online) Dependence of the integrated num- 
ber of eigenvalues N\ with |A| > 0.25 (red/gray squares) and 
A I > 0.1 (black circles) as a function of the total number 
of processes TV for versions of Linux Kernels. The values of 
iV correspond (in increasing order) to Linux Kernel versions 
1.0, 1.1, 1.2, 1.3, 2.0, 2.1, 2.2, 2.3, 2.4 and 2.6. The power 
law iV A oc TV" has fitted values ^| A |>o.25 = 0.622 ± 0.010 and 
^|a|>o.i — 0.630 ±0.015. Inset shows data for the Google ma- 
trix G* with inverse link directions, the corresponding expo- 



nents are v* x[>0 _ 2S 



0.696±0.010 and v\ 



|A|>0.1 



0.652±0.007. 



The eigenvalues \ and the right eigenvectors 'ipi(j) sat- 



isfy the equation V ■, Gjj>ipi(j') = Xiipi(j). To obtain 
the spectrum of the PCN with large sizes N we used 
the Arnoldi method from ARPACK library [1J] that al- 
lowed to find all eigenvalues with |A| > 0.1 up to the 
maximum N = 285509. The complex spectrum of A; 
is shown in Fig. Q] for three versions of Linux Kernel. 
With the increase of the version number and the size N 
certain characteristic structures appear in the distribu- 
tions of Xi = \Xi\ exp(iipi) in the complex plane. Thus, 
there are clearly visible lines at real axis and polar an- 
gles ip = 7r/2,27r/3,47r/3,37r/2. The later are related 
to certain cycles in procedure calls, e.g. an eigenstate at 
X{ = 0.85exp(i27r/3) is located only on 6 nodes. In Fig.Q] 
we also present the spectrum of the Google matrix G* 
which is obtained by an inversion of link directions in the 
adjacency matrix of the PCN as it was proposed in (23J . 
It is characterized by a similar spectrum structure but a 
larger radial distribution of A; . 

The dependence of the normalized integrated density 
of states W(7), with eigenvalues < -fi < 7, on the 
eigenvector relaxation rate 7 is shown in Fig.[2j Here 7$ = 
—2 In |Aj I and W{"f) is normalized on the number of states 
N\ in the interval A < |Aj| < 1 with A = 0.1. In average 
the dependence W(-f) is the same for all 10 versions of 
Linux Kernel even if there are certain fluctuations. There 
are clear vertical steps in W(j) at -f m = — 21n(a/m) 
which are due to a presence of many degenerate states at 
the corresponding values of A (e.g. 116 for A = a for the 
version 2.6.32). Such degeneracy is known to be typical 
of the WWW networks (see e.g. 
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FIG. 4: (Color online) Cluster growing fractal dimension of 
the Linux Kernel directed networks. The average mass (M c ), 
with uniformly distributed seed, is shown for different Kernel 
versions as a function of the network distance I, in logarithmic 
scale. The corresponding fractal dimensions d for all versions 
are shown in the inset panel, where the power law fit is done 
for the range 1 < I < 10. The dashed and dot-dashed lines (in 
both panels) represent the power law dependence (M c ) oc l d , 
corresponding to fractal dimensions d — 2v — 1.30 and d = 
2v* — 1.26 from the fractal Weyl law of Fig. [3] respectively. 



The network size N grows with the version number 



of Linux Kernel corresponding to its evolution in time. 
The dependence of N\ on N, shown in Fig. [31 clearly 
demonstrates the validity of the fractal Weyl law with 
the exponent v « 0.63 for G (we find v* « 0.65 for G*). 
We take the values of v for A = 0.1 where the number of 
eigenvalues N\ gives a better statistics. Withing statisti- 
cal errors the value of v is not sensitive to the cutoff value 
at small A. The matrix G* has slightly higher values of v. 
These results show that the PCN of Linux Kernel has a 
fractal dimension d — 2v w 1.26 for G and d = 1v pa 1.3 
for G*. 



III. FRACTAL DIMENSION OF LINUX 
KERNEL NETWORK 

To check that the fractal dimension of the PCN in- 
deed has this value we compute the dimension of the 
network by another direct method known as the cluster 
growing method (a description of this method can be find 
in |2l[). In this method the average mass or number of 
nodes (M c ) is computed as a function of the network dis- 
tance I counted from an initial seed node with further 
averaging over all seed nodes. In a dimension d the mass 
(M c ) should grow as (M c ) oc l d that allows to determine 
the value of d for a given network. 

It is important to note that up to now the cluster grow- 
ing method (see e.g. |2l|) was used only for undirected 
networks. Our network is a directed one and, since we 
did not find relevant methods for such a case in a liter- 
ature, we generalized this method in a simple way using 
the network distance I following only outgoing links. The 
average of (M c (l)} is done over all nodes. Due to global 
averaging the method gives the same result for the matrix 
with inverted link direction (indeed, the total number of 
outgoing links is equal to the number of ingoing links). 
We will see that the fractal dimension obtained by this 
generalized method is very different from the case of con- 
verted undirected network. 

The dependence of (M c ) on I, obtained for the directed 
network, is shown in Fig. 01 Up to a saturation regime 
related to a finite number of nodes the growth of (M c ) 
is satisfactory described by a power law with a certain 
fractal dimension d. The straight lines, with the value 
of d obtained from the fractal Weyl law, give in average 
a satisfactory description of the growth (M c (l)). More 
detailed data, with the fitted values of d for the interval 
1 < I < 10, are shown in the inset panel of Fig. @] for var- 
ious network sizes N. For Linux versions up to V2.3 the 
fractal dimension remains in the interval 1.3 < d < 1.45 
with an average d « 1.4. However, we have d ps 1.6 for 
V2.4 and d w 1.9 for V2.6 versions. We attribute the 
strong deviation for the version V2.6 to the well known 
fact that significant rearrangements in the Linux Kernel 
have been done after version V2.4 [25, 26]. Also the data 
of Fig. 2] clearly show the saturation of (M c (l)) at rela- 
tively small values (M c (l)) < 200 compared to the whole 
size of the network N < 285509. This happens due to 
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FIG. 5: (Color online) Cluster growing function for the undi- 
rected version of the Linux Kernel networks. The average 
mass (M c ), with uniformly distributed seed, is shown for dif- 
ferent undirected Kernel versions as a function of the network 
distance I, in logarithmic scale. 
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FIG. 6: (Color online) Localization properties of Google ma- 
trix eigenstates for Linux Kernel versions 2.0.40 2.4.37.6 and 
2.6.32 with N = 14079, 85756, 285509 respectively. The par- 
ticipation ratio £ of eigenstates is shown as a function of mod- 
ulus of the corresponding eigenvalues |A|; the values of < £ >, 
averaged over all given interval of |A|, are 9.5, 8.6, 17.3 respec- 
tively. 



existence of many isolated communities in the network 
which become linked only due to regularization of the 
Google matrix ([IJ . This is the reason of strong degener- 
acy of the eigenvalue X = a which has 116 values for the 
version V2.6. Also, in contrast to undirected networks, in 
directed networks there are nodes with only ingoing links 
that give an important reduction of the average mass. 

Thus in view of the above restrictions we consider that 
there is a rather good agreement of the fractal dimen- 
sion obtained from the fractal Weyl law with d ~ 1.3 
and the value obtained with the cluster growing method 



which gives an average d ~ 1.4. The fact that d is ap- 
proximately the same for all versions up to V2.4 means 
that the Linux Kernel is characterized by a self-similar 
fractal growth in time. The closeness of d to unity signi- 
fies that procedure calls are almost linearly ordered that 
corresponds to a good code organization. Of course, the 
fractal Weyl law gives the dimension d obtained during 
time evolution of the network. This dimension is not nec- 
essary the same as for a given version of the network of 
fixed size. However, one can expect that the growth goes 
in a self-similar way [lj| [22| and that static dimension 
is close to the dimension value emerging during the time 
evolution. This can be viewed as a some kind of ergodic- 
ity conjecture. Our data show that this conjecture works 
with a good accuracy up to the Linux Kernel V.2.6. 

Here it is important to note that the properties of the 
directed network are very different from the converted 
undirected network with the same nodes. The later is 
obtained from the directed network by replacing all di- 
rected links by undirected ones. The results obtained by 
the clustering growing method for this undirected net- 
work are shown in Fig. \5[ We find that it is difficult to 
fit the dependence (M c (l)) by a simple power law in this 
case. A very rough estimate gives very high value d ~ 5 
being by a factor 4 larger than the dimension of the di- 
rected network. Also the saturation level of (M c ) is now 
comparable with the whole system size iV while for the 
directed network it is smaller than N by three orders of 
magnitude. 

Usually the methods of computation of fractal dimen- 
sion have been developed for undirected networks (see 
e.g. H3, EHUmi). Our example shows that the di- 
rected networks may have rather different characteristics 
compared to their undirected versions. Thus new meth- 
ods for computation of fractal dimension in directed net- 
works should be developed. Our results show that the 
generalized version of the method described in [21| work 
well in the case of directed networks. 



IV. GOOGLE MATRIX EIGENSTATES 

Before we discussed the properties of the spectrum of 
the Google matrix of the PCN. Let us now analyze the 
properties of its eigenstates tpi(j). As it is usual for the 
disordered systems with the Anderson localization (see 
e.g. [2!|), we use the PArticipation Ration (PAR) £ to 
characterize the localization properties of the eigenstates. 
It is defined as & = (^ l^(j)| 2 )7£, \Mj)\ 4 - Physi- 
cally, £j gives the number of nodes effectively occupied by 
an eigenstate tpi(j). The dependence of £ on |A| for three 
versions of the Linux Kernel is shown in Fig. [51 There is 
only a small increase of £ with size N which is however 
much smaller than the increase of N from one version to 
another. Also even maximum values of £ remain by a 
factor 5 x 10 3 smaller than the system size N. Thus we 
conclude that the eigenstates are essentially localized. 

The dependence of probability of eigenstates, with 




FIG. 7: (Color online) Coarse-grained probability distribution 
\ipi(j))\ 2 for the eigenstates of the Google matrix for Linux 
Kernel version 2.6.32. The horizontal lines show the first 64 
eigenvectors ordered vertically by |Aj| from |A| = 1 on bottom 
to |A| ~ 0.4 on top (only one state is shown for degenerate 
eigenvalues, all states are normalized by £V \ipi(j)\ 2 = 1, in- 
dex j corresponds to the order given by the PageRank located 
at the bottom line). The top panel shows the coarse-grained 
complete network of 285509 sites with 307 cells, each con- 
taining 930 sites; the bottom panel shows the coarse-grained 
distribution for first 300 cells, each containing 62 sites (bar 
and arrow on the top panel mark the region of bottom panel) . 
Probability is proportional to color with black for zero and 
yellow/white for maximum. 



largest values of |Aj|, on the PageRank index j is shown 
in Fig. [7]for the version 2.6.32. The main peaks of prob- 
ability are scattered over the whole system size N. The 
appearance of vertical stripes is clearly visible. This sig- 
nifies that statistically the positions of peaks are repeated 
from one eigenstate to another. It is possible that this 
effect is related with certain cycles between procedures 
in the Linux Kernel. These cycles, or communities as it 
is used to say for the WWW, appear even on a far tail 
of the PageRank. The origin of such weakly isolated cy- 



cles and correlations between different relaxation modes 
requires further detailed studies. 



V. DISCUSSION 

Here we characterized the directed network by the 
Google matrix and investigated its spectrum and eigen- 
state properties. We think that such an approach gives 
much richer charctcristization of networks. Indeed, as it 
was shown in [20], the spectrum of the Google matrix 
depends in a sensible way on network structure so that 
networks with rather similar distribution of links may 
have very drastic differences in the spectrum (e.g. pres- 
ence or absence of gap in a vicinity of A = 1 ) . We think 
that further studies of the Goggle matrix properties of 
complex networks will bring us deeper understanding of 
their generic properties. 

In conclusion, our results show that the Google matrix 
of the PCN of the Linux Kernel is characterized by the 
fractal Weyl law with the fractal Weyl exponent v ps 0.65. 
This value corresponds to the fractal dimension of the 
network d « 1.3. This dimension is smaller than two 
and thus the fractal Weyl law becomes well visible. In 
contrast to that for networks with the dimension d > 2 
the fractal Weyl law is replaced by a usual dependence 
N X (xN since v = d/2>l. It is known that the WWW 
has the fractal dimension d « 4.1 [2lJ that explains why 
the fractal Weyl law does not work for the WWW univer- 
sity networks discussed in [2fJ] . On the basis of our result 
we argue that the fractal Weyl law should appear in all 
directed networks with the fractal dimension d < 2. It is 
important to note that the fractal Weyl exponent v is not 
sensitive to the exponent characterizing the decay of the 
PageRank: indeed, according to [23| the later remains 
the same for the WWW and the PCN of Linux Kernel 
while the values of v are different. We think that the 
further studies of the fractal Weyl law for directed net- 
works with d < 2 will bring us to a deeper understanding 
of their scale-free properties. Our results also show that 
the properties of directed networks can be rather differ- 
ent from their undirected versions. 
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